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基于 径 问 基 小 波 网 络 的 二 代 棉 铃 虫 
卵 峰 日 预测 模型 
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摘要 : 为 了 提高 害虫 预报 的 准确 率 , 将 径 回 基 小 波 网 络 首次 引入 农作物 害虫 预测 预报 领域 , 改进 了 径 向 基 小 波 网 络 
的 学 习 算 法 , 使 之 适合 于 害虫 预测 的 应 用 : 利用 径 向 基 小 波 函 数 族 时 、 频 两 域 支撑 完全 或 部 分 覆盖 被 分 析 数 据 序列 
时 、 频 两 域 支撑 的 原理 来 确定 小 波 函 数 族 尺度 参数 和 平移 参数 取 值 ; 根据 中 心 向 量 之 间 的 欧式 距离 大 小 来 初步 第 选 
隐 含 层 神经 元 。 在 实例 分 析 中 , 本 文 利 用 1966 - 1995 年 山东 省 惠 民 县 棉铃 虫 Helicoverpa armigera 的 监测 数据 建立 了 
基于 径 向 基 小 波 网 络 的 2 代 棉 铃 忠 孵 量 峰值 日 期 预测 模型 利用 1996 - 2000 年 的 监测 数据 对 模型 进行 了 检验 。 检 
验 结 果 表 明 : 在 5 年 的 预测 数据 中 , 4 年 的 预测 数据 偏差 在 3 d 以 内 , 男 外 1 年 的 预测 数据 偏差 4d, 预测 效果 令 人 满 
意 。 本 文 为 害虫 预测 预报 研究 提供 了 一 种 可 行 的 新 方法 。 
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Forecasting model for the oviposition peak day in the second generation of 
Helicoverpa armigera ( Lepidoptera: Noctuidae ) based on radial basis 


Wavelet network 
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University, Nanjing 210095, China; 2. Plant Protection Station of Shandong Province, Jinan 250100, China) 
Abstract : To improve the accuracy of crop pest forecasting, this paper introduced and applied radial basis 











wavelet network into the area of crop pest forecasting for the first time. The author modified the learning 
algorithms of radial basis wavelet network for application in pest forecasting. The scale and translation 
parameters were determined by the theory that time-frequency support of analyzed data sequence is 
covered with time-frequency support of radial basis wavelet functions. Based on the Euclidean distance 
between central vectors, the hidden-layer neurons are selected preliminarily. At case study, the 
investigation data of Helicoverpa armigera in Huimin, Shandong between 1966 and 1993 were used to 
establish the forecasting model of oviposition peak day in the second generation of H. armigera based on 
radial basis wavelet network, while the investigation data between 1996 and 2000 were used to test the 
model. The test results showed that the forecasting deviation of four years was less than three days and the 
forecasting deviation of one year was four days. The forecasting results proved satisfactory. This paper 
developed a new studying method for crop pest forecasting. 
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如 何 提 高 害虫 预测 的 准确 率 是 农作物 病虫害 预 
测 预 报 工作 中 的 核心 和 热点 问题 之 一 。 田 间 害 虫 种 
群 的 消长 和 影响 害虫 生长 发 育 的 外 界 因素 共同 组 成 
了 一 个 复杂 的 非 线性 系统 , 因此 , 采用 非 线 性 预测 
分 析 方 法 进行 害虫 的 预测 预报 具有 传统 线性 预测 方 
法 无 法 比拟 的 优点 。 小 波 神 经 网 络 (wavelet neural 


network, WNN) (Zhang and Benveniste, 1992; Baskshi 
and Stephanopoulous, 1993; Pat and Krishnaprasad, 
1993 ) 是 融合 了 小 波 分 析 和 神经 网 络 2 种 非 线 性 理论 
的 数学 建 模 方法 , 具有 很 强 的 非 线 性 逼近 能 力 。 但 是 ， 
WNN 在 高 维 输入 的 情况 下 会 产生 “ 维 数 灾 难 ”, 这 是 
因为 小 波 基 萄 数 会 随 着 维 数 的 增加 而 呈 指 数 级 的 增 
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长 , 导致 网 络 结构 和 运算 代价 急剧 增加 。 为 此 , 很 多 
学 者 (Zhang et al.,，1995 ; 吴 钦 军 等 ,1996; Zhang， 
1997; 万 建 等 ,2002 ) 分 别 从 不 同 的 角度 探讨 了 如 
何 有 效 减 少 高 维 WNN 的 神经 元 数量 以 及 降低 运算 
代价 , 取得 了 一 定 的 效果 , 但 并 未 从 根本 上 避免 “ 维 
数 灾难 ”。Kugarajah 和 Zhang (1995 ) 构造 了 多 维 单 
尺度 径 问 小 波 框 染 , 给 有 效 解决 高 维 WNN 所 面临 的 
问题 开辟 了 一 条 新 的 途径 。 由 于 径 向 小 波 对 维 数 不 
敏感 , 因此 , 以 多 维 单 太 度 径 问 小 波 框架 为 基础 的 径 
问 基 小 波 网 络 便 可 以 有 效 避 开 “ 维 数 灾难 ”问题 。 

目前 , 多维 径 回 基 小 波 网 络 在 国内 应 用 较 少 ， 
关 清 等 (2003 ) 将 其 应 用 到 脑 电 逆 问 题 求解 ,而 在 害 
忠 预 测 领 域 尚 无 应 用 先例 。 害 忠 预 测 模 型 通常 具有 
多 维 或 高 维 输入 的 特点 , 适合 于 建立 多 维 径 向 基 小 
波 网 络 , 基于 此 , 本文 将 多 维 单 尺度 径 癌 基 小 波 网 
络 引 入 害虫 预测 领域 , 对 算法 进行 了 适当 的 调整 和 
修改 , 使 之 适合 于 害虫 预测 方面 的 应 用 , 并 通过 对 
1996 - 2000 年 山东 惠 民 县 2 代 棉 铃 虫 卵 峰 日 的 预 
测 来 检验 预测 模型 的 效果 。 
1 建 模 方法 
1.1 多 维 径 向 基 小 波 网 络 结构 

径 癌 基 神 经 网 络 (radial basis function neural 
network ，RBFNN ) 是 采用 递减 径 回 基 函 数 做 为 隐 合 
层 激励 函数 的 前 馈 神 经 网 络 , 一 个 d 维 输入 1 维 输 
出 的 RBFNN 结构 如 图 1 所 示 。 

网 络 输出 为 : 

f(x) = 之 WiD( lx -cs |) + 

其 中 , nn 为 隐 含 层 神经 元 的 数量 , W 表示 隐 含 
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图 1 径 回 基 神 经 网 络 的 结构 
Fig. 1 The structure of radial basis neural network 

X,: 神经 网 络 第 i 个 输入 变量 The ;th input variable of neural network ; 
C;: 隐 含 层 第 ;个 神经 元 激励 函数 的 中 心 向 量 Central vector of 
stimulation function of the ith neuron in hidden layer; 中: 作为 隐 含 层 激励 
函数 的 径 向 基 羡 数 Radial basis function as stimulation function of hidden 
layer; | ，|| : 欧 几 里 得 范 数 Euclidean norm; 玩 : 隐 含 层 第 i 个 神经 元 
与 输出 层 之 间 的 连接 权 值 Weight between the ith neuron of hidden layer 
and output layer; 6: 输出 层 偏 移 量 Bias of output layer; 3: 求 和 孔 数 
Summation function; y: 神经 网 络 输出 Output of neural network. 


层 第 个 神经 元 与 输出 层 的 连接 权 值 , @(，) :BR 
一 R 作为 隐 含 层 激励 函数 的 径 向 基 范 数 ,c; sR 表 
未 隐 含 层 第 个 神经 元 激励 函数 的 中 心 同 量 , 5 为 
输出 层 偶 移 量 。 

用 多 维 单 广 度 径 癌 小 波 基 郴 数 族 取 代 RBFNN 
隐 含 层 的 激励 函数 就 可 以 构造 出 多 维 单 广度 径 癌 基 
小 波 网 络 , 伸缩 步 长 为 2、 位移 步 长 为 1 的 d 维 输 
入 工 维 输出 的 径 同 基 小 波 网 络 结构 如 网 2 所 示 。 
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图 2 多 维 单 尺度 径 向 基 小 波 网 络 结构 
Fig. 2 The structure of radial basis wavelet network based on multi-dimension and single-scaling 
X,: 神经 网 络 第 i 个 输入 变量 The ;th input variable of neural network; m,: 尺度 参数 Scale parameter; n,: 位 移 向 量 Translation parameter; 中 : 径 


问 小 波 函 数 Radial wavelet function; | 。 


: 欧 几 里 得 范 数 Euclidean norm; WW;: 隐 含 层 第 i 个 神经 元 与 输出 层 之 间 的 连接 权 值 Weight between 


the ;th neuron of hidden layer and output layer; 6b: 输出 层 偏 移 量 Bias of output layer; 3: 求 和 函数 Summation function; y: 径 向 基 小 波 网 络 输出 


Output of radial basis wavelet network. 
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网 络 输出 为 : 


f(x) = DAC - nn) + bb = 


mal len 1),, 
i 2 


式 中 , 于 为 隐 含 层 第 字 个 神经 元 结 点 与 输出 层 
的 连接 权 值 , m; 为 尺度 参数 , n; 为 位 移 问 量 , 中 为 
径 问 小 波 函 数 , 下 为 径 向 小 波 孔 数 的 多 维 形 式 , 6 
为 输出 层 偏 移 量 ,s 为 隐 含 层 神 经 元 结 点 数量 。 

本 文中 的 径 癌 小 波光 数 采用 文献 (Kugarajah 
and Zhang，1995 ) 中 提 到 的 墨西哥 草帽 函数 
( Mexcian Hat, Mexh ) 。 
1.2 建 模 步骤 

普通 RBFNN 的 训练 过 程 分 为 2 步 : 第 1 步 为 
无 教师 学 习 , 确定 训练 输入 层 与 隐 含 层 间 的 权 值 ; 
第 2 步 为 有 教师 学 习 , 确定 训练 隐 合 层 与 输出 层 间 
的 权 值 和 侦 移 量 。 径 癌 基 小 波 网 络 的 训练 过 程 与 普 
通 RBFNN 基本 相同 , 但 是 算法 不 一 样 。 吴 清 等 
(2003 ) 详细 探讨 了 径 回 基 小 波 网 络 的 学 习 算 法 , 本 
文 作者 对 其 中 不 适合 于 害虫 预测 应 用 的 具体 实现 过 
程 进行 了 改进 , 建立 了 径 癌 基 小 波 网 络 害 中 预测 模 
型 的 训练 算法 , 包括 以 下 3 步 : 

第 1 步 (Step 1) : 确定 尺度 参数 m 和 位 移 参 数 
n 的 取 值 范围 , 初步 得 到 隐 含 层 神 经 元 的 构成 。 

从 图 2 可 以 看 出 , 要 确定 隐 含 层 神 经 元 、 求 出 
输入 层 与 隐 含 层 连接 权 值 ,必须 要 先 确 定 尺度 参数 
和 位 移 问 量 的 取 值 , 吴 清 等 (2003 ) 根 据 妹 奎 斯 特 采 
样 定理 来 确定 尺度 参数 的 取 值 范围 ,而 对 于 害虫 预 
测 应 用 来 说 , 调查 数据 本 喘 就 是 离散 数据 ,所 以 无 


法 通过 样本 数据 的 采样 间隔 来 确定 尺度 参数 。 本 文 
作者 对 其 进行 了 改进 , 提出 采用 时 域 和 频 域 交叉 、 
履 善 的 方法 ,具体 来 说 ， 就 是 在 建立 预测 模型 时 ， 
作为 隐 含 层 激 励 函 数 的 小 波 基 函数 的 时 域 和 频 域 支 
撑 区 间 至 少 能 够 履 关 到 1 个 样本 上 后， 从 而 对 未 知 映 
射 的 估计 过 程 产 生 贡 献 。 改 进 算法 中 的 关键 点 是 分 
析 确 定 所 有 输入 因子 的 周期 频率 变化 范围 , 可 采用 
带宽 参数 为 1、 中 心 频率 为 1 的 Morlet 复 小 波 分 析 
确定 , 分 析 方 法 参考 王 文 圣 等 (2003 ) 的 方法 。 

第 2 步 (Step 2) : 在 每 一 个 上 广度 空间 中 挑选 能 
较 好 地 对 样本 进行 解释 的 小 波 基 函数 ,进一步 往 选 
第 1 步 得 到 的 隐 含 层 神 经 元 。 

第 3 步 (Step 3): 用 施 密 特 逐步 正 交 化 方法 去 
掉 见 余 的 隐 含 层 神经 元 , 求解 最 佳 小 波 基 于 集 ， 同 
时 完成 网 络 训练 。 

本 研究 全 部 数据 运算 采用 Matlab 数值 分 析 软 
件 编 程 完 成 。 


2 ”实例 分 析 


2.1 实例 资料 和 数据 处 理 

实例 模型 为 利用 气象 因素 预测 2 代 棉 铃 虫 卵 量 
峰值 日 期 (简称 卵 峰 日 ) 。2 代 棉 铃 虫 虫 情 资 料 来 源 
于 山东 省 农作物 病 虫 预测 预报 观测 数据 集 , 采用 的 
是 山东 惠 民 县 1966 - 2000 年 2 代 棉 铃 虫 发 生 期 实 
况 汇总 数据 ; 气象 资料 来 源 于 山东 省 气象 局 1966 - 
2000 年 逐日 降水 资料 、 逐日 气温 资料 、 逐 旬 平 均 气 
温 资料 、 逐 月 平均 气温 资料 、 逐 旬 平 均 降 水 资料 、 
逐 月 平均 降水 资料 。 其 中 , 1966 -1995 年 30 年 的 


表 1 山东 惠 民 县 2 代 棉 铃 虫 卵 峰 日 径 向 基 小 波 网 络 模型 的 输入 变量 和 输出 变量 


Table 1 The inputs and outputs for radial basis wavelet network model for the oviposition peak day 


in the second generation of Helicoverpa armigera in Huimin, Shandong 


代号 Symbol 意义 Meaning 
11 4 月 中 旬 平 均 气 温 Average temperature of middle April 
LD 4 月 下 旬 平 均 气 温 Average temperature of late April 
B 5 月 平均 气温 Average temperature of May 
14 6 月 上 旬 平 均 气 温 Average temperature of early June 
输入 变量 、 
5 6 月 中 旬 平 均 气 温 Average temperature of middle June 
Input variables 16 4 月 中 旬 降 水 量 Rainfall of middle April 
17 4 月 下 旬 降 水 量 Rainfall of late April 
I8 5 月 降水 量 Rainfall of May 
I9 6 月 上 旬 降 水 量 Rainfall of early June 


输出 变量 


| Output 
Output variable 


2 代 棉 铃 虫 卵 峰 日 


Oviposition peak day of the second generation of Helicoverpa armigera 


模型 输出 变量 为 2 代 槐 铃 虫 卵 量 峰 值 日 期 , 把 5 月 31 日 作为 基准 值 0, 其 他 日 期 按 步 长 值 1 依次 增 减 。The 


output variable of model is the oviposition peak day of the second generation of H. armigera. May 31 was taken as the zero 


base value and the other date changed with each day as step 1. 
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实况 资料 作为 模型 训练 样本 ,1996 -2000 年 5 年 的 
实况 资料 作为 模型 测试 样本 , 用 来 检验 模型 效果 。 
本 研究 选择 了 9 个 气象 因子 作为 模型 的 输入 ， 
模型 输出 为 2 代 卵 峰 日 (以 5 月 31 日 为 0, 其 他 日 
期 依次 类 推 )( 表 1)。 
为 了 保证 良好 的 计算 性 能 ， 降低 运算 的 复杂 
度 , 对 样本 数据 做 归 一 化 处 理 , 处 理 原则 是 : 


original - location 
data = 一 一 一 一 一 一 — 
scale 


其 中 ，original 为 原始 数据 ，location 为 每 一 变 
量 样本 数据 中 的 最 小 值 ，scale 为 样本 数据 的 最 大 值 
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减 去 最 小 值 。 
2.2 数据 分 析 和 建 模 

首先 , 根据 时 域 和 频 域 尾 盖 的 原则 来 确定 径 问 
基 小 波 函 数 伸缩 参数 和 位 移 向 量 的 取 值 。 模 型 输入 
因子 的 样本 数据 经 过 归 一 化 处 理 后 , 所 有 因子 的 样 
本 数据 都 在 0 到 1 之 间 变 动 , 因此 , 所 有 输入 因子 
的 时 域 文 撑 范 围 都 是 [0,1] 。 

利用 Morlet 复 小 波 进行 小 波 变 换 , 分 别 计算 9 
个 输入 因子 的 小 波 方差 , 并 绘制 各 因子 的 小 波 尺度 
方差 图 ( 王 文 圣 等 , 2003 ) 。 
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图 3 输入 变量 小 波 方差 图 
Fig. 3 Wavelet variance figure of input variables 
11 -9 所 代表 意义 同 表 1。The meaning of 11 -Jo is the same as in Table 1. 


通过 分 析 各 因子 小 波 尺 度 方 差 图 , 可 以 得 到 9 
个 目 变 量 输入 因子 的 周期 范围 和 相应 的 频率 变动 范 
围 , 取 各 目 变 量 输入 因子 频率 变化 范围 的 并 集 , 得 

27T 2 万 
到 | 和 5 ， 3 | 

Mexh 阴 数 的 频 域 有 效 文 撑 范 围 可 根据 频率 域 


能 量 集中 的 原理 ( 张 邦 礼 等 ,1995 ) 来 确定 ,本 研究 
取 占 总 能 量 99. 9% 的 区 域 作为 母 小 波 函 数 有 效 频 
域 支 撑 范 围 , 为 [10,3] 。 径 向 基 神 经 网 络 中 的 径 向 
基 哨 数 必须 是 正定 的 (Buhmann，1998), 径 问 基 小 
波 网 络 中 的 Mexh 小 波 函 数 也 必须 满足 这 一 条 件 ， 
因此 , 本 文中 的 Mexh 母 小 波 函 数 时 域 支 撑 范 围 确 
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定 为 | -1,1]。 伸 缩 步 长 为 2、 尺 度 参数 为 m、 位 移 
步 长 为 1、 位移 参数 为 n 的 Mexh 小 波 函 数 族 的 时 
频 文 撑 可 以 表示 为 : 

OO =L2"( -1+n),2"(1+n)] x[0,3 x2"| 

根据 Mexh 小 波 函 数 的 频 域 支撑 需要 与 输入 因 
子 频率 变化 范围 相互 覆盖 或 互相 交叉 的 原理 , 则 3 
x2 ”>27/[2S , 可 求 得 m <3.577,， 因为 m 为 整 
数 , 所 以 m 最 大 为 3。 而 由 于 Mexh 小 波 函 数 正 频 
域 最 小 为 0, 则 实际 上 尺度 参数 m 可 以 取 为 任意 比 
3 小 的 整数 (包括 0 和 负 整 数 ) 。 但 是 , 对 于 尺度 参 
数 m 范围 的 确定 , 需要 以 最 小 的 范围 能 够 覆盖 因子 
频率 范围 为 原则 ,所 以 对 于 尺度 参数 最 小 值 的 确定 
采用 以 下 计算 方法 : 3 x 2”=27w/3, 则 m= 
-0.518, 因此 尺 度 参 数 m 最 小 为 0， 尺 度 参数 mm 
的 取 值 范围 确定 为 [0,3], 即 尺度 参数 可 取 值 为 : 
0,1,2,3,。, 

当 尺 度 参 数 m 为 0 时 , 为 了 使 Mexh 小 波 函 数 
族 时 域 支 撑 能 够 完全 禾 盖 输入 因子 时 域 支 撑 范 围 ， 

-1l+n>0,1+n<l 

此 , 尺度 参数 m 为 0 时, 位 移 参数 取 值 范围 
为 [0,1], 由 于 位 移 步 进 为 1, 则 位 移 参数 取信 为 0 
和 1。 由 于 各 维 输入 因子 都 经 过 了 归 一 化 处 理 , 变 
动 范围 一 致 , 所 以 , 各 维 位 移 参 数 取 值 相同 , 都 为 
[0,1]。 同 理 , 可 以 求 得 尺度 参数 为 1, 2, 3 时 ,位 
移 参 数 取 值 也 为 [0,1 |]。 
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建 模 第 2 步 是 隐 含 层 神 经 元 的 初步 第 选 。 根 据 
上 一 步 所 求 得 的 位 移 参 数 的 取 值 , 不 同 输入 维 的 位 
移 量 可 以 组 成 不 同 尺度 下 隐 含 层 神经 元 激励 函数 的 
中 心 问 量 , 每 一 个 中 心 右 量 都 对 应 着 相应 的 隐 含 层 
神经 元 , 这 样 获得 的 初始 神经 元 数量 在 高 维 的 情况 
下 会 比较 庞大 , 因此 , 我 们 需要 在 保证 能 够 完全 和 覆 
六 数 据 信号 时 频 空间 的 基础 上 精简 基 函 数 的 数量 ， 
减少 基 消 数 的 元 余 性 。 

同一 尺度 参数 下 ,两 个 基 泡 数 的 中 心 向 量 之 间 
欧式 距离 如 果 过 小 , 我 们 认为 2 个 基 函 数 相似 ， 则 
它们 从 样本 数据 中 提取 的 信息 也 十 分 相似 , 可 以 只 
取 其 一 。 本 例 中 , 经 过 第 1 步 运 算 , 伸缩 参数 为 0， 
1,2 和 3 时 , 分 别 可 以 组 合 出 512 个 中 心 向 量 。 根 
据 备 选 向 量 的 多 少 , 本 例 设 定 当 2 个 中 心 疝 量 之 间 
的 欧式 距离 小 于 0. 34 时 , 则 表示 由 两 中 心 同 量 组 成 
的 基 函 数 相 似 , 可 以 排除 挥 1 个 基 消 数 。 经 过 沛 
选 , 从 2 048 个 基 函 数 中 初步 第 选 出 1 024 个 向 量 。 

最 后 , 将 30 个 训练 样本 输入 , 采用 上 自 适 应 正 交 
投影 算法 进行 网 络 学 习 训 练 , 最终 将 隐 含 层 神 经 元 
数量 确定 为 30 个 , 并 计算 出 隐 含 层 与 输出 层 连 接 
权 值 以 及 输出 层 偏 移 量 。 

至 此 , 我 们 最 终 构造 了 一 个 具有 30 个 隐 含 层 
神经 元 方太 、9 维 输 入 、1 维 输出 的 径 回 基 小 波 网 
络 , 利用 1996 -2000 年 5 年 的 实际 监测 数据 对 网 
络 模型 进行 检验 , 表 2 给 出 了 检验 结果 , 模型 预报 


表 2 模型 输出 与 发 生 实况 比较 


Table 2 Comparison between the outputs of model and the real results 


标准 化 输入 因子 Standardized input factors 


年 份 

Year 11 了 B 14 BS 16 
1996 0. 153 0.727 0.445 0.513 0.495 0 
1997 0. S01 0. 525 0.134 0.029 0.930 0 
1998 1.000 0.463 0.221 0 0.607 0 
1999 0.689 0.581 0.253 0.789 -0.518 0. 193 
2000 0.614 0.568 0.292 0.440 1.000 0.133 





转化 为 日 期 
输出 资料 
Date of 
17 18 I9 Output eal 
forecasting 
date 
0.097 0.030 ”0.006 0. 560 6/26 6/26 
0.210 0.482 0.002 -0.075 6/19 6/21 
0.313 0.967 0.221 —0.082 6/19 6/20 
0.078 0. 545 0. 065 0.714 6/28 6/24 
0.011 0.134 0.004 -0.083 6/19 6/22 


I ~ 19 所 代表 意义 同 表 1, 所 有 变量 均 经 过 标准 化 处 理 , 模型 输出 同样 为 标准 化 值 , 模型 输出 值 再 经 过 反 标 准 化 处 理 并 加 上 基准 值 
5 月 31 日 后 得 到 预测 日 期 , 反 标 准 化 公式 : Forecasting date = output x (31 -20) +20. The meaning ofIl -JI9 is the same as in Table 1. All 


input variables were processed by standard treatment and the output of model was standard value too. Date of forecasting was gotten as a result of 


the sum of output anti-standardized and the base value May 31. Anti-standardized formula: forecasting date = output x (31 - 20 ) +20. 
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与 讨论 


径 回 基 小 波 网 络 作为 一 种 融合 WNN 和 RBFNN 
的 神经 网 络 , 兼 具 WNN 通 近 能 力 强 (Zhang and 
Benveniste ，1992 ; 李 力 和 方 华 京 , 2000) 和 RBFNN 
训练 速度 快 、 不 存在 局 部 最 优 ( 周 俊 武 等 , 2001; 共 
小 红 等 ,2006 ) 的 优点 , 同时 , 克服 了 WNN 在 高 维 
输入 时 所 产生 的 “ 维 数 灾难 ”问题 , 大 大 简化 了 高 维 
WNN 的 算法 。 对 本 例 而 言 , 奉 采 用 张 量 积 的 形式 
构建 多 维 小 波 基 ， 则 在 达到 同样 的 精度 下 ， 建 立 
WNN 所 需要 的 权 值 数目 是 以 指数 倍 N (CN 是 指 对 
一 维 函 数 盘 近 时 所 需要 的 权 值 数 目 ) 的 规模 来 衡量 
的 , 即使 N 是 个 很 小 的 整数 , 所 需要 的 权 值 数 目 也 
将 很 庞大 , 这 也 是 WNN 多 用 于 一 维 或 低 维 函 数 通 
近 的 原因 ， 而 本 文采 用 的 径 癌 基 小 波 网 络 在 9 维 输 
入 的 情况 下 最 终 只 需要 30 个 隐 含 层 神 经 元 即 可 完 
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成 建 模 , 极 大 地 缩减 了 模型 的 规模 。 

农作物 病虫害 的 发 生 和 和 危害 作为 一 种 生态 系 
统 , 具有 复杂 的 非 线 性 特点 ， 因此 常 用 的 统计 预测 
在 预测 精度 上 难以 达到 满意 的 效果 ,小波 网 络 结合 
了 小 波 变换 和 神经 网 络 2 种 非 线 性 理论 , 在 模拟 非 
线性 系统 上 具有 一 定 的 优势 。 表 3 是 基于 误差 反问 
传播 的 神经 网 络 (back-propagation neural network ， 
BPNN) 和 RBFNN 对 本 例 的 预测 结果 , 2 种 神经 网 
络 都 采用 3 层 结 构 ，BP 网 络 隐 含 层 结构 分 别 设 定 
为 18, 19, 20 和 30 个 神经 元 ，RBFNN 的 spread 参 
数 取 值 为 1.5。 将 表 3 与 表 2 对 比 后 发 现 , 针对 同 
样 数据 进行 建 模 , 径 向 基 小 波 网 络 的 预测 精度 要 高 
于 BPNN 和 RBFNN。 在 建立 的 4 种 BPNN 中 , 预测 
精度 最 高 的 是 隐 含 层 具 有 19 个 神经 元 的 BPNN ,而 
RBFNN 和 径 回 基 小 波 网 络 一 样 都 具有 30 个 隐 含 层 
神经 元 ,这 意味 着 径 癌 基 小 小 网络 继承 了 普通 
RBFNN 所 具有 的 结构 比 BPNN 略 复杂 的 特点 。 


表 3 BPNN 和 RBFNN 预测 结果 
Table 3 Forecasting of BP network and RBF network 


BPNN 
Wd 18 神经 元 19 神经 元 20 神经 元 30 神经 元 RBFNN 站 
Eighteen neurons Nineteen neurons Twenty neurons Thirty neurons 
1996 6/26 6/24 6/22 6/22 6/24 6/26 
1997 6/23 6/22 6/24 6/28 6/22 6/21 
1998 6/22 6/22 6/26 6/28 6/23 6/20 
1999 6/27 6/24 6/26 6/23 6/29 6/24 
2000 6/26 6/25 6/25 6/26 6/19 6/22 


BPNN: 基于 误差 反 向 传播 的 神经 网 络 Back-propagation neural network; RBFNN: 径 向 基 神 经 网 络 Radial basis function neural network ; 
表 中 日 期 (月 /日 )( 除 最 后 一 栏 ) 为 每 种 模型 对 1996 -2000 年 山东 惠 民 县 2 代 棉 铃 虫 卵 峰 日 预测 结果 。The dates (month/day) (except 


the last column ) in the table are the forecasting results of oviposition peak day of the second generation of Helicoverpa armigera between 1996 and 


2000 im Huimin, Shandong, China. 


不 论 是 统计 预报 模型 , 还 是 神经 网 络 等 非 线 性 
预报 模型 ,都 是 在 春 干 年 发 生 实况 监测 数据 的 基础 
上 进行 总 结 、 归 纳 ， 从 而 得 出 发 生 特 点 ， 进 而 建立 
模型 ， 因 此 , 模型 的 建立 要 求 历 史 数 据 越 多 越 好 。 
历史 实况 数据 越 多 , 模型 越 能 够 总 结 出 病虫害 发 生 
的 真实 特点 ， 从 而 模型 的 准确 率 越 高 。 

病虫害 的 发 生 受 多 种 因素 影响 , 包括 气象 条 
件 、 病 虫 基数 、 作 物 长 势 和 天 敌 等 因素 , 因此 ， 如 
何 选取 模型 的 输入 因子 是 建立 模型 时 要 面 对 的 首要 
问题 。 生 态 系统 内 部 各 因素 间 具 有 复杂 的 相互 关 
系 , 所 以 , 仅仅 依据 相关 系数 或 侦 相 关系 数 等 统计 


学 上 的 相关 理论 来 判断 和 选择 输入 因子 是 狭 隆 的 相 
关 理 论 , 而 非 线 性 相关 理论 在 目前 尚未 有 较 合理 的 
计算 方法 , 因此 如 何 确定 模型 输入 目 变 量 和 输出 因 
变量 之 间 的 非 线性 相关 性 也 是 以 后 需要 着 重 研 究 的 
问题 。 本 文 在 研究 过 程 中 , 也 曾 尝 试 过 利用 小 波 分 
解 得 到 的 小 波 系 数 来 计算 两 因素 间 的 相关 性 , 以 后 
可 以 做 深入 的 人 研究。 
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